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The non-equilibrium dynamics of an isolated quantum system after a sudden quench to a dy¬ 
namical critical point is expected to be characterized by scaling and universal exponents due to the 
absence of time scales. We explore these features for a quench of the parameters of a Hamiltonian 
with 0{N) symmetry, starting from a ground state in the disordered phase. In the limit of infinite 
N, the exponents and scaling forms of the relevant two-time correlation functions can be calculated 
exactly. Our analytical predictions are confirmed by the numerical solution of the corresponding 
equations. Moreover, we find that the same scaling functions, yet with different exponents, also 
describe the coarsening dynamics for quenches below the dynamical critical point. 
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I. INTRODUCTION 

Understanding the non-equilibrium dynamics of iso¬ 
lated macroscopic quantum systems is currently one of 
the most intriguing challenges of statistical physics [T]. 
The interest in this problem has been recently revived by 
impressive advances in the physics of ultracold atoms 12 - 
[^, which paved the way to the experimental investiga¬ 
tion of the real-time dynamics of almost isolated quan¬ 
tum many-body systems iHISl with tunable properties 
such as interactions, dimensionality, internal symmetries, 
etc. Consequently, what was an academic question until 
very recently, i.e., whether and how an isolated quantum 
system thermalizes after a quench of the parameters of 
its Hamiltonian, became a subject of extensive investi¬ 
gation [IMH]. Remarkably, the study of the approach 
to thermal equilibrium disclosed a number of novel fea¬ 
tures such as the existence of prethermal quasi-stationary 
states, predicted in Ref. cni, and subsequently experi¬ 
mentally observed in Refs. HOI HD- It was then ar¬ 
gued |22II!??T| that this prethermalization occurs due to 
the system retaining memory of the integrable part of 
the Hamiltonian, which, in the absence of non-integrable 
terms, would lead the system towards a non-thermal 
state. The question of whether these prethermal states 
can be represented by a generalized Gibbs ensemble is 
also a topic of current, active research |31It55] . 

It has also been argued that a quantum system can un¬ 
dergo a dynamical phase transition (DPT) |3HH46] in the 
prethermal regime, with some of the features of an equi¬ 
librium quantum phase transition, such as the appear¬ 
ance of diverging length- and time-scales. In particular, 
such a DPT was shown DU mi to be related, in some 
respects, to an equilibrium one with an effective temper¬ 
ature proportional to the amount of energy injected by 
the quench into the system. On this basis, one expects 
some degree of universality to emerge close to a DPT, 
with spatial and temporal correlations characterized by 
critical exponents which do not depend on the micro¬ 


scopic details of the system. In particular, the lack of 
time scales combined with the breaking of translational 
symmetry in time due to the sudden quench is expected 
to cause an algebraic relaxation of the system to its sta¬ 
tionary state, a phenomenon which has been referred to 
as quantum aging, in analogy with the slow dynamics ob¬ 
served in glasses [15]. While aging is known to occur in 
classical critical systems USED], only recently the same 
issue has been investigated after quenches in open |51b 
ESI or isolated m quantum systems. However, while in 
the former case the critical point responsible for aging 
is the thermal one dictated by the presence of a ther¬ 
mal bath, in the latter the critical point has an intrinsic 
non-equilibrium nature. 

The consequences of this difference has been recently 
explored in Ref. m, where a system with vector order 
parameter cj) with N components and a Hamiltonian with 
0{N) symmetry was studied by means of a perturbative 
renormalization-group approach. Aging was shown to oc¬ 
cur in an early stage of the evolution, when the system 
exhibits prethermalization, rather than in the eventual 
thermalized state. This prethermal state was recognized 
to be similar to the non-equilibrium stationary state ap¬ 
proached by the same model in the limit A —>■ oo, in 
which it becomes integrable |42fHl] . This fact suggests 
that the 0(A)-symmetric Hamiltonian for A —>■ oo is 
the integrable Hamiltonian responsible for the prether¬ 
malization occurring for finite A. 

In the present work we report a complementary non- 
perturbative analysis of the problem analyzed perturba- 
tively in Ref. m , by considering the post-quench dynam¬ 
ics in the limit A —> oo, when the model is numerically 
and analytically exactly solvable. 

The structure of the paper is as follows: in Sec. 
we introduce the model and the protocol of the quench, 
and derive the relevant self-consistent evolution equa¬ 
tions which are exact in the limit A —>■ cx). In Sec. nil 
we discuss their analytic solution in the absence of time- 
and length-scales in the system, i.e., after a quench at 
or below the dynamical critical point. In particular, we 
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provide analytical predictions for the relevant scaling ex¬ 
ponents and scaling functions for generic spatial dimen¬ 
sionality d of the model. In Sec. |IV| we present the results 
of the numerical solution of the evolution equations for 
a quench to the critical point (Sec. IV AI and below it 
(Sec. |I VB ), finding excellent agreement with the ana¬ 
lytic predictions; for comparison we briefly discuss the 
latter case in the context of coarsening (Sec. IVCl. Fi¬ 
nally, in Sec. [V] we summarize our findings and briefly 
discuss their possible implications for recent experimen¬ 
tal proposals. 


II. MODEL AND QUENCH PROTOCOL 

We consider a system described by the following 0{N)- 
symmetric Hamiltonian in d spatial dimensions 


the quartic interaction in Eq. (Q decouples as |421 - HH IM] 

^ 2{N + 2)((/)2)</,2 _ N{N + (2) 

as shown in Ref. [^, this decoupling corresponds to the 
Hartree-Fock approximation, which becomes exact for 
fV —)■ 00 . Once inserted into Eq. Q, the dynamics of 
the various components of the fields decouples and each 
of them is ruled (up to an inconsequential additive con¬ 
stant) by the effective time-dependent quadratic Hamil¬ 
tonian 

Hes{t) = ^ / [n^ + + reff(t)02] , (3) 


where rgg(t) is determined by the condition 

7'eff(t) = r-k j^((/>2(x,t)). (4) 


H{r, u) 




u 




( 1 ) 

where (p = {pi ,..., (/)jv) is a bosonic held with N compo¬ 
nents, while n is the canonically conjugated momentum 
with [(()a(x), nb(x')] = — x')6ab- Note that only 

the scalar products pP = (p ■ (p = 11^ enter 

Eq. 0 . as a consequence of the symmetry. The coupling 
M > 0 controls the strength of the anharmonic interac¬ 
tion. At time t < 0 the system is prepared in the disor¬ 
dered ground state \po) of the pre-quench non-interacting 
Hamiltonian H(Qg,0), while the system is evolved with 
the post-quench Hamiltonian H{r, u) for t > 0. This pro¬ 
tocol was studied in Ref. m while in Refs. |4()l[i2Vlii] u 
had the same non-zero value in the pre- and post-quench 
Hamiltonian. For the latter, it was shown that the sys¬ 
tem undergoes a dynamical phase transition (DPT) at a 
certain r = Tc, below which coarsening occurs. As ar¬ 
gued in Ref. m, a non-vanishing u in the present initial 
Hamiltonian would only renormalize the value Qq and 
therefore we expect the model discussed here to undergo 
a DPT such as the one studied in Refs. |40l l42lj44] . The 
upper critical dimension of the DPT is dc = 4, as in 
the corresponding equilibrium theory at finite tempera¬ 
ture muz]. Accordingly, the critical, collective proper¬ 
ties for d > 4 are completely described by the correspond¬ 
ing non-interacting theory, as discussed in Sec. |HI[ For d 
smaller than the lower critical dimensionality di = 2 |44| . 
instead, this DPT is not expected to occur. However, 
since a counterpart of the Mermin-Wagner theorem is 
not known out of equilibrium, the occurrence of a differ¬ 
ent DPT cannot be ruled out for d < di. 

In the following we will exploit the fact that, if the 
0{N) symmetry in the initial state is not broken, the 
average (papb) (where (...) = (^o| • ■ • \'Po)) vanishes un¬ 
less a = b, with its non-vanishing value independent of a, 
and equal to the fluctuation {p'^) of a generic component 
p of the held. In the limit N —b oo, the model can be 
solved by taking into account that, at the leading order. 


In a field-theoretical language, reff(t) plays the role of a 
renormalized square effective mass of the held p. Due to 
the quadratic nature of ideff, it is convenient to decom¬ 
pose the held p into its Fourier components p^, according 
to p{x,t) = / d'^/c ())k(t)e®'^'^/(27r)'^, with an analogous 
decomposition for H. Each of these components can be 
written in terms of the annihilation and creation oper¬ 
ators HSJ SU Ok and aj^, respectively, diagonalizing the 
initial Hamiltonian: 

P\,{t) = fk{t)ai, + /k(t)aLk, (5) 

where /k(t) = /_k(t) is a complex function depending on 
both time t and momentum k. Note that the canonical 
commutation relations between p and H imply |43| 

2Im[/k(t)/^(t)] = l. (6) 

The Heisenberg equation of motion for p\^ derived from 
the Hamiltonian 0 yields the following evolution equa¬ 
tion for /k(t): 


/k + + res{t)]h = 0, (7) 

with k = |k|. This equation is supplemented by the initial 
conditions 

/k(0) = l/y/2uJok and /k(0) = -iVwofc/2, (8) 

where Wq^. = + ^Iq, which can be obtained by diag¬ 

onalizing the quadratic pre-quench Hamiltonian and by 
imposing the continuity Eauzi of 0k(i) at < = 0. 

Since the Hamiltonian is quadratic, all the information 
on the dynamics is encoded in its two-time functions, 
such as the retarded (Gn) and Keldysh (Gk) Green’s 
functions, which are defined by |58| 

6u.-k'iGR{k,t,t') = 9{t - t'){[pi^{t),pu’{t')]), (9) 

Sk-k>iGK{k,t,t') = {{pk{t),pkt{t')}), (10) 
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where 5k.k' = (27r)‘^5^‘^^(k — k'), while 9{t < 0) = 0 and 
9{t > 0) = 1. Accordingly, reff in Eq. Q can be expressed 
in terms of Gk as 


11 f d ^ A* 

res(t) = r+— J j^^iGK(k,t,t)h(klK), (11) 

where the function 


h{x) 


0 for cc 1, 
1 for X 1, 


( 12 ) 


implements a large-k cut-off at the scale A in order to 
make the theory well-defined at short distances, while it 
does not affect it for fc <C A. By using Eq. (|^ in Eqs. ( [I0| ) 
and ([^, we find that Gk,r can be written in terms of the 
function /k as: 


iGK{k,t,t')=2Re[Ut)fl{t')], (13) 

GR{k,t,t') = 29{t (14) 


The dynamics of the system can be determined by solv¬ 
ing the set of self-consistency equations given by Eqs. ([^, 
(111 and (13). Generically, these equations do not admit 
an analytic solution and therefore one has to resort to 
numerical integration. Nevertheless, in the following sec¬ 
tion, we show that some quantities can be analytically 
calculated when the post-quench Hamiltonian is close to 
the dynamical critical point. 


III. DYNAMICAL PHASE TRANSITION AND 
SCALING EQUATIONS 

In Refs. |42lH4] it was shown that, after the quench, 
the system approaches a stationary state, in which the 
effective Hamiltonian becomes time-independent. Such 
a stationary state was then argued to be non-thermal 
as a consequence of the integrability of the model. In 
particular, in Ref. B], it was demonstrated that when 
r in Eq. 0 is tuned to a critical value the long¬ 
time limit r* of the corresponding effective parameter 
reff(t) in Eq. 0 vanishes and therefore the fluctuations 
of the order parameter become critical, signalling the oc¬ 
currence of a dynamical phase transition. More precisely, 
as a consequence of the divergence of the spatial corre¬ 
lation length ^ = (r*)“^/^, the correlation functions at 
long times acquire scaling forms characterized by univer¬ 
sal critical exponents. 

Similarly to the case of classical systems after a quench 
of the temperature EHI, the correlation functions exhibit 
dynamical scaling forms not only in the steady state, but 
also while approaching it m-- in particular, relying on 
dimensional analysis and on the lack of additional time- 
and length-scales at r = one expects the effective 
value Tesit) to scale as: 

ree{t) = ^cr(At), (15) 


where the function a is normalized by requiring a(oo) = 
1, such that reff(t) vanishes at long times as 


reg(t) = ^ for At » 1, 


(16) 


while a is a dimensionless quantity. The non-universal 
corrections introduced by a (At) — 1 to this long-time 
limit are negligible for At ^ 1. On the contrary, for 
At 1, they become dominant and non-universal be¬ 
havior is displayed. Accordingly, one can identify a mi¬ 
croscopic time m tA ~ A ^ which separates these two 
regimes: for 0 < t < t^ the dynamics is dominated by 
non-universal microscopic details; for t ^ t\, instead, the 
dynamics becomes universal. This discussion assumes 
that the function cr(r) has a well-defined limit as r —>■ oo, 
which might not be the case in the presence of oscillatory 
terms. In fact, as shown in the numerical analysis pre¬ 
sented in Sec. |IV| the non-universal function a depends 
on how the cut-off A is implemented in the model, i.e.. 


on the choice of the function h(x) in Eq. (111. In par¬ 


ticular, the choice of a sharp cut-off turns out to make 
ct(t) oscillate, masking the universal long-time behavior 

reff(t) ~ 


As a consequence of the universal form of Eq. (16) for 


t > tA) the correlation functions are expected to exhibit 
scaling properties. In order to show this, it is convenient 
to rescale time and write the function /k(t) as /k(t) = 
gk(kt). Inserting Eq. (16) into Eq. one finds the 

equation for g-k{x): 


9kix) + (l + ;^) 9k(x) = 0, 
valid for x = kt > kt^, whose solution is: 


(17) 


9k{x) = Vx[AkJa(x) -I- BkJ_a(x)], (18) 

where Jq(x) is the Bessel function of the first kind and 


a = 


— a. 


(19) 


Below we show that it is consistent to assume a < 1/4 
and therefore a to be real. The constants Ak and Hk in 


Eq. (181 are fixed by the initial condition of the evolution. 


as discussed below. For later reference, we recall that 


j (2;') ^ J + “)’ a:<l, 

|cos(a: — a7r/2 — 7 r/ 4 )-^/ 27 ( 7 ray, a; ^ 1, 

( 20 ) 

where r(a:) is the Euler gamma function [59]. Note that 
Eq. ( jlSl encodes the complete dependence of /k(t) on 
time t for t ^ t\, whereas its dependence on the wave 
vector k is encoded in the yet unknown functions Ak and 
Hk- By using the following identity for the Wronskian of 
Bessel functions 


Ux)j'_^(x) - J.Ux)J'^(x) = 


( 21 ) 
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one can show that Eq. ^ requires the coefficients 
Bk to satisfy the relation: 


Im[^kS^] = - 


4sin(a7r) k 


( 22 ) 


This relation is not sufficient in order to determine com¬ 
pletely Ak and Sk unless the full functional form of reff(t) 
is taken into account, including its non-universal behav¬ 
ior for t < ^a; this would allow us to fix Ak and i?k on the 
basis of the initial conditions for the evolution, which at 


present cannot be reached from Eq. (181, it being valid 


only for t ^ However, for a deep quench — such 
as that one investigated in Ref. m — with fig 
the initial conditions ([^ for the evolution of /k become 
essentially independent of k and read: 

/k(0) ~ 1/v^, /k(0) -^VHo72. (23) 

At time t ~ fA ~ A“^, /k(tA) can be calculated from a se¬ 
ries expansion = fk{0)+tAfk{0)+t\jiki0)+O{t\) 

and by using Eqs. ([^ and ( |23[ ) one may readily conclude 
that its dependence on k comes about via (fc/A)^ and 
A:/(figA); accordingly, at the leading order, it can be ne¬ 
glected for k A. On the other hand, /k at can be 
evaluated from Eq. (18 1 and, since kt\ ~ k/A ^ 1, we 
can use the asymptotic form of the Bessel functions for 
small arguments, finding: 


/k(tA) OC Ak(fcfA)^^^^“ + Bk{ktAy^^~ 


(24) 


In order to have /k(tA) independent of k at the leading 
order, it is then necessary that: 

Ak ~ A(fc/A)-i/2-“, Rk ^ B(A:/A)-i/2+“, (25) 

where A and B are yet unknown complex numbers. 
While this scaling is expected to be true for k A, 
non-universal corrections may appear for A: ~ A. Note 
that Eq. (25) is consistent with (22), provided that: 


lm[AB*] = - 


ttA ^ 

4sin(Q;7r) 


(26) 


The numerical analysis discussed in Sec. |IV| actually 
shows that Ak turns out to be purely imaginary, while 
Rk real. Combining Eqs. (25) with (18), (13), and (14), 


one finds a simple form for the retarded Green’s function 


GR{k,t,t') = -d{t-t') 




2sina7r ' (27) 

[Ja{kt)J-aikt') — J-a{kt)Ja{kt')] , 

while the one for iGK{k,tA') is somewhat lengthy and 
thus we do not report it here explicitly. In order to deter¬ 
mine a and therefore a from the self-consistent condition 


in Eq. (11), it is actually sufficient to know iGxikAA') 
for t = t', which is given by 

iGxik, t, t) = 2At { |A|2(A:/A)-2“ J2(A;A) 

-k |R|2(A:/A)2“ J2^(fci) + 2Re[AB*]J^{kt)J-c,{kt)} . 

(28) 


Accordingly, while Eq. (22) is sufficient in order to deter¬ 


mine the complete form of Gr, the one of Gk still con¬ 
tains unknown coefficients A and B which are eventually 
determined by the initial conditions; nevertheless, the 
scaling properties of both of these functions are already 
apparent. In fact, their dynamics is characterized by two 
temporal regimes, which we refer to as short (kt <C 1) 
and long (kt 3> 1) times. Stated differently, the tem¬ 
poral evolution of each mode k has a typical time scale 
~ k~^ determined by the value of the momentum it¬ 
self and the corresponding short-time regime extends to 
macroscopically long times for vanishing momenta. 

Let us focus on Gr in Eq. (27): at short times t' < 
t it becomes independent of k 




1/2-c 




(29) 


where we used the asymptotic expansion in Eq. (20). For 


well-separated times ti$>f, the second terms in brackets 
is negligible and GR{k,t S> t',t') displays an algebraic 
dependence on the ratio t'/t. At long times k~^ ^ t',t, 
instead, Gr becomes time-translational invariant and by 
keeping the leading order of the asymptotic expansion of 
the Bessel functions in Eq. (|20|), it reads 


GRik,t,t')^-0it-t' f’'^l (30) 

which is nothing but the Gr of a critical Gaussian Hamil¬ 
tonian after a deep quench m- Similarly, at short times 
t <C A;“^ and to leading order in At, Gxik, t, t) reads (see 


Eqs. (28) and (20)) 


Ol —2a I A |2 

zGx(A:,t,t):^ I ' (At)^+^- 
1 "'(I -I- a) 


(31) 


which is independent of the momentum k, as it is the 
case for Gr within the same temporal regime. At long 
times t ^ k~^, instead, one finds 

^GKik,t,t)^^(^^'^^ , (32) 

to leading order in k/A, where the possible oscillating 
terms have been neglected, as they are supposed to av¬ 
erage to zero when an integration over momenta is per¬ 
formed. The resulting expression turns out to be time- 
independent and, contrary to what happens with Gr, it 
does not correspond to the Gk of a critical Gaussian the¬ 
ory after a deep quench m, reported further below in 
Eq. (33). This fact should be regarded as a consequence 


of the non-thermal nature of the stationary state which is 
eventually reached by the system and which retains mem¬ 
ory of the initial state. Since the effective Hamiltonian 
(§ is Gaussian, it is tempting to interpret the anomalous 
momentum dependence ~ k~^~^°‘ in Eq. (32) in terms of 


a quench in a truly Gaussian theory. In the latter case, 
it was shown that, for deep quenches, 

iGK{k,t,t) ~ 


(33) 
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which is similar to an equilibrium distribution with an ef¬ 


fective temperature Tes ~ flo- Accordingly, Eq. (321 can 


be regarded as resulting from a quench of a Gaussian the¬ 
ory with a momentum-dependent initial “temperature” 
noik) ~ 


Equations (311 and (32) show that the term propor¬ 
tional to |Ap in the expression (281 of Gfcik, t,t) is dom¬ 


inant at both long and short times: accordingly, we can 
neglect the remaining ones and the Keldysh Green’s func¬ 
tion with two different times t, t' acquires the scaling 
form 


iGK{k,t,t') = 2|A|^ 


1+2q; 


Vk'^tt' Ja{kt)Ja{kt'), 


(34) 

which can be derived from Eqs. (IlSL (25), and (13). This 


expression can now be used in order to derive from the 
self-consistency equation (111 the value of the constant 
a and therefore of the exponent a (see Eq. @) which 
characterizes the scaling forms ( [2^ and (3^. For At ^ 1 
one can approximate refi(t) with Eq. (161 and, by assum¬ 
ing a sharp cut-off 


h{x) = 0(1 — x), 


Eq. (Ill becomes 


a 

t2 


ad 


\A\'^uA‘^RdA^t) 


(35) 


(36) 


at the critical value Tc of r. In this expression we intro¬ 
duced the function 

Rd,a{x)=^f (37) 

-'O 

while Qd = fld/{‘^T^Y with 12^ = (T{dj2) the d- 

dimensional solid angle. Expanding for large ar¬ 

gument a; 3> 1, one finds that 


Rd,a(.x) = Wd,a{x) + + 


Ji) 


.(2) 


d-l-2a 


O 


(38) 

where Wd,a(x) is a fast oscillating function which is a 
consequence of the sharp cut-off considered in the in¬ 
tegral over momenta and are certain coefficients, 
the relevant values of which are provided further below. 


Once the expansion (381 is plugged into Eq. (36), the 


r.h.s. of the latter can be expanded in decreasing powers 
of At which have to match the term on its l.h.s., resulting 
in a set of conditions fixing the values of Tc, a, and u. 
Notice that the oscillations contained in Wd^a{x) would 
be compensated by sub-leading terms in the l.h.s. of 
Eq. (|^ (which are not reported, see the discussion af¬ 
ter Eq. @) , thus confirming the non-universal nature of 
Wd^a{x)- The value of is determined such as to cancel 
the constant contribution ex on the r.h.s., while a has 
to be fixed such that = 0 in order to cancel the term 

a^ot 


cx (At)^+2“-d which cannot be matched by the l.h.s. of 
Eq. ( [M| |61] . This procedure — typically used for solv¬ 
ing this kind of self-consistency equations |491 15^ 155] — 
can be regarded as a systematic way of canceling terms 
which result from corrections to scaling and, therefore, it 
allows a comparison with the results obtained within a 
renormalization group approach (with some exceptions, 
see Ref. [53 )• For our purposes, it is sufficient to focus 
on the condition 


,(i) 

‘-'d,a 


1 




2v^r(l-f+a)r(l-f+ 2 a) 


= 0, (39) 


which is solved by requiring the argument of one of the 
two r functions in the denominator to equal a non¬ 
positive integer. The two corresponding infinite sets of 
solutions S' and S for a are given by 


c' j d 2 
D = < —z- n 


S = 


n>0 


d — 2 n 


, (40) 


ra >0 


with integer n. The physically relevant solution can be 
selected by requiring a to match its Gaussian value a = 0 
(or, equivalently, a = 1/2, see Eq. at the upper crit¬ 
ical dimensionality dj. = 4 of the model UH133 Note that 
this Gaussian value of a can be easily inferred by inspect¬ 
ing the scaling behavior of Gr^k for the quench towards 
a non-interacting theory (compare, e.g., Eqs. (321 and 
(331). Accordingly, one finds a single possible solution 
from S, i.e.. 


d-2 


a = 


and 


d ( d 
4 I 4 


(41) 


with a > 0 for all values of d between di =2 and dc = 4, 
while for d > 4 the Gaussian theory applies and therefore 


1 


a = 


and 


a = 0. 


(42) 


( 2 ) 

For d > dc, oc a vanishes and the leading temporal 


dependence of the r.h.s. of Eq. (361 at long times is oc 
4 ^^^i/ 2 (At)“('^“^) and therefore 


reff(t) oc (At) 


-{d-2) 


(43) 


up to oscillating terms, instead of the behavior oc t ^ 
in Eq. (16). Equations (41) and (42) together with 


Eqs. (27) and (34) completely characterize the scaling 


behavior of Gk,r after a deep quench to the critical 
point. These results can be compared with the predic¬ 
tions of Ref. m, formulated with a dimensional expan¬ 
sion around the upper critical dimensionality dc of the 
0{N) model, which is expected to reduce to the present 
one for N —)■ oo. There, the time dependence of Gk,r 
for kt, kt' 1 and t' t was parametrized in terms of 
the exponent 0 as 


Gfl(fc, t, t') oc —t(t'/t)® and Gk oc {tt') 


i-e 


(44) 
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from which it follows that 9 is related to the exponent 
a introduced in Eqs. (|2^ and p4| by 0 = 1/2 — a. 
Moreover, it was found in Ref. |47| that 9 = e/4 + ©(e^) 
with e = 4 — d, in agreement with Eq. (411, which yields 
6» = 1 - d/4 for d < 4. 

We note here that among the remaining solutions of 
Eq. (39) in iS and S' which do not match the Gaussian 


value at d = dc, only one turns out to be compatible 
with having a positive value of the coupling constant u 
in Eq. (36). This solution belongs to S', and is given by 


d- 2 


<^co — 


. , (3 — d)(d — 1) , 

with flco = - 1 -; (45) 


remarkably, it turns out to be related to the coarsening 
occurring after a quench to r < rc, as we argue and 
demonstrate further below. 

A comparison between these results and those for the 
corresponding exactly solvable classical quench is dis¬ 
cussed in Appendix [A| 


IV. NUMERICAL RESULTS 


In order to test the quality of the analytical predictions 
of the previous section, we studied in detail the numeri¬ 
cal solution of the evolution equations for /k(t), un¬ 
der the constraint provided by Eq. with Gk given 

by Eq. (13). The numerical integration of these equa¬ 


tions has been performed using an algorithm based on 
the Bulirsch-Stoer method |^, while the integrals over 
momentum k have been computed using the extended 
Simpson’s rule |M] with a mesh of 7.5 x 10^ points. 

In Sec. |IV A| we consider the case of a quench to the 
critical point, comparing the numerical results for the 
relevant correlation functions with the analytical predic¬ 
tions derived in Sec. |III[ In Sec. |IV B| we present, in¬ 
stead, results for a quench below the critical point, and 
show some numerical evidence of the emergence of scaling 
properties during coarsening. 


A. Quench to the critical point 


In order to determine the critical value of the pa¬ 
rameter r in Eq. (IR one can conveniently use the ansatz 
proposed in Ref. gllMI, i e.. 


rc = 


u 

4! 


dfik 2k^ + Qg 


h{klA), 


(46) 


which turns out to correctly predict Tc also beyond the 
case of a deep quench. The rationale behind this ansatz 
relies on the assumption that Gi<-(fc,t, t) at t —>■ oo has 
approximately the same form as for the case of a quench 
to r = 0 in the non-interacting case m = 0 [SJ. Although 
the validity of this assumption for a quench to the critical 
point is questionable because, as shown in Sec. |III[ Gk 
differs significantly from the non-interacting case u = 0 


(see, e.g.. Ref. SZl), Eq. ( |46[ ) anyhow provides accurate 
predictions for the value of re- 

The accuracy of the ansatz (161 for the long-time be¬ 
havior of reff(t) can be tested by calculating reg(t) ac¬ 
cording to Eqs. (Ill and (13), based on the numerical 
solution of the evolution equation 0 for /k(t). 

In Fig. we show the time dependence of reff(t) after a 
deep quench. In particular, the upper left panel demon¬ 
strates, in spatial dimension d = 3.4, that while reg(t) at 
long times approaches a finite value for r > rc, it generi- 
cally vanishes for both r = rc and r < In particular, 
the corresponding decay turns out to be ~ in both 
cases, as indicated by the dashed lines, in agreement with 
the ansatz (161 at criticality. Such an algebraic behavior 
actually sets in after some time ~ C1(A“^), where A 
is the cut-off employed in the algorithm. As expected, 
the actual possibility of detecting this algebraic decay 
depends on the way the model is regularized, i.e., on the 
specific function h{x) used in Eq. (Ill in order to in¬ 
troduce the cut-off A. This is illustrated by the upper 
central panel of Fig. [^for a quench at criticality r = Tc 
in d = 3, in which the momentum integral in Eq. (Ill 
is regularized with a sharp cut-off as in Eq. ( |35[ ) (upper¬ 
most curve) or with a smooth exponential h[x) = e”’’’ 
or Gaussian h{x) = e““ A functions (lower curves), with 
characteristic scale A. (Note that the value of is also 
affected by this choice, according to Eq. (461.) While 
the sharp cut-off causes persistent oscillations in |refi(t)|, 
which mask the expected behavior ^ the smooth 
ones are qualitatively similar and they reveal this alge¬ 
braic decay after some time ^a. The persistent oscilla¬ 
tions displayed for a sharp cut-off function are expected 
to be sub-leading compared to the A-independent decay 
~ in a formal expansion in decreasing powers of A: 
accordingly, their amplitude is expected to decrease as 
A increases. This is clearly demonstrated by the curves 
in the upper right panel of Fig. which, from top to 
bottom, slowly approach the expected algebraic behav¬ 
ior (dashed line) upon increasing the value of A. Accord¬ 
ingly, in order to detect the universal behavior ~ 
in the presence of a sharp cut-off, very large values of A 
have to be used, resulting in a longer computational time 
with respect to that required by an exponential or Gaus¬ 
sian cut-off function h. The onset of a scaling regime 
for reff(t) for a quench at criticality is also expected to 
be influenced by the val ue o f the pre-quench parameter 
as discussed in Sec. |III| and in Ref. m- In partic¬ 
ular, while the analytic investigation in Sec. |III| assumes 
a deep quench, i.e., Qq 3> A, it is instructive to check 
numerically how the actual value of Qq influences the 
time Ia after which the expected universal algebraic be¬ 
havior ~ sets in. The curves in the lower left panel 
of Fig. [2 show, from top to bottom, that increases 
significantly upon decreasing the value of Oq, until the 
eventual algebraic behavior is completely masked by the 
initial non-universal transient occurring at t < <a when 
Qq ^ A. If the spatial dimensionality d of the model 
is larger than the upper critical dimensionality dc = 4, 
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the leading-order temporal decay of refi(t) is no longer 
proportional to because the corresponding propor¬ 
tionality constant a vanishes (see Eq. In this case, 

reff(t) still vanishes at long times with the algebraic law 
^ ^-(ii- 2 ) given in Eq. (43l. This dependence is shown 


in the lower central panel for d = 4.1, 4.3 and 4.5. The 
theoretical prediction for the exponent of this decay is 
indicated by the corresponding dashed lines, while its 
prefactor is a non-universal constant which depends — 
in contrast to the case d < 4 — on the actual values of 
the parameters of the system, e.g., flo, as shown in the 
lower right panel for d = 4.5. 


As discussed in Sec. Ill the scaling behavior of reff(t) 


implies the emergence of algebraic dependences on time 
in Gji Ki numerical evidences for which are presented in 
Fig. 15 for a deep quench occurring at d = 3. In partic¬ 
ular, GK{k,t,t) at criticality is expected to display the 
short-time scaling in Eq. ( |M] ) for t k~^ (which actually 
extends to long times for fc = 0), while in the long-time 
limit t ^ k~^, Gk displays (up to oscillating terms) the 
scaling in Eq. (321 as a function of fc. In order to corrob¬ 


orate these predictions, the left and right upper panels 
of Fig. [5 show the dependence on time t and momen¬ 
tum k, respectively, of Gk within these two regimes, for 
quenches occurring at r > Tc (lowermost curves), r = rc 
(intermediate curves), and r < (uppermost curves). 
The upper left panel demonstrates that Gxik = 
after a quench to criticality, grows in time as ~ (in¬ 
termediate dashed line) in agreement with Eqs. and 
for d = 3. If the quench occurs above criticality, 
instead, Gk displays oscillations with an asymptotic pe¬ 
riod oc ^ = (r*)“^/^ where r* = reff(t —>■ oo) as expected 
on general grounds m and as suggested by the low¬ 
ermost curve. Remarkably, an algebraic behavior ^ 
(uppermost dashed line) emerges for quenches occurring 
below the critical point (i.e., with r < rc), which signals 
that the corresponding coarsening occurs on a dynamical 
length scale which grows in time [BS]. Further below (see 
Fig. m, we discuss the dependence of the power of this 
algebraic decay on the spatial dimensionality d. We an¬ 
ticipate here that the present law ^ agrees with what 
one would obtain by extending the scaling prediction in 
Eq. (311 to quenches below Tc and by using the value 

instead of a. The upper right 


Q^co = 1/2 from Eq. 
panel, instead, reports the dependence on k of the value 
eventually reached by GK^k, t, t) at a fixed but long time 
t^k~^. In particular, at criticality r = (intermediate 
curve), Gk approaches, up to oscillatory terms, the alge¬ 
braic behavior ~ k~^/'^ (lower dashed line) in agreement 
with Eqs. (321 and ( [41] ) in d = 3. For r < rc (uppermost 
curve), Gk still displays, up to oscillatory terms, an alge¬ 
braic dependence on fc, but with a different power ~ k~^ 
(upper dashed line) which again agrees with the extension 
of the critical scaling form (32) below rc with a replaced 
by ctco- When the quench occurs, instead, above the crit¬ 
ical point, Gk tends to a constant, up to oscillations. In 
all the cases illustrated in Fig. [5 non-universal contribu¬ 
tions affect the various curves for fc > A = 7r/2, due to the 


effects of the regularizing function h. As far as (fc, t, t') 
is concerned, Eq. (27) provides its complete expression 


within the scaling regime at criticality r = re- In par¬ 
ticular, Gfl(fc = 0,t,t') acquires the scaling form (29), 
while GK{k,t,t') becomes time-translationally invariant 
at long times k~^ t' < t as in Eq. ( |M| . The lower left 
panel of Fig. [^ shows that Gji{k = 0,t,t')/t with t' < t 
becomes indeed a function of the ratio t'/t only. At criti¬ 
cality (lowermost curve) this agrees with what is expected 
on the basis of Eqs. (p^) and (411 in d = 3 (which renders 
a = 1/4), with |Gi^(fc = 0,t,r)/t| oc $i/ 4 (F/f) (dashed 
line), where we define 


$a(a:) = X 


l/2-a _ 2,1/2 +o 


(47) 


For a quench below rc, instead, the same quantity be¬ 
comes |Gi^(fc = 0,t,t')/t| (X 1 — {f /t) (dashed line); this 
agrees with /t), i.e., with what one would infer by 

extending the critical scaling function (291 below Cc and 
by using the value aco = 1/2 in Eq. (45) for the expo¬ 


nent a. The inset shows, instead, the numerical data for 
Gnik = after a quench to r > rc, as a function 

of t' and fixed t > t', which is characterized by persis¬ 
tent oscillations. As mentioned above, for this kind of 
quench, reff(t) approaches a finite asymptotic value r* at 
large times and, up to the leading order, the system be¬ 
haves as a Gaussian model quenched at r*, for which the 


response function is given by [57] Eq. (30) with 


\/k^ 


+ r* 


(48) 


The dashed line reported in the inset of the figure, which 
is actually indistinguishable from the numerical data, 
corresponds to this theoretical prediction with k = 0 
(see also next panel), confirming its accuracy. The lower 
right panel of Fig. |5 shows Gfi{k,t > t',f) as a func¬ 
tion of k{t — t') for two fixed long times t and t' < t 
and upon varying k. The main plot shows the corre¬ 
sponding numerical curves both for r = rc and r < rc, 
which are however perfectly superimposed and practi¬ 
cally indistinguishable from the theoretical prediction in 
Eq. (30) (dashed line), the latter being independent of 


the actual value of a. The inset, instead, shows Gb. for 
r > rc'- also in this case, the numerical data are indistin¬ 
guishable from the corresponding theoretical prediction 
obtained on the basis of Eqs. (30) and (48), as explained 


above while illustrating the previous panel. 

The numerical results presented in Figs. |5 and [prefer 
to a quench of the model in a certain spatial dimension¬ 
ality d. In order to test both the predictions in Eq. 
and some of the features of the scaling functions at criti¬ 
cality we repeated the analysis of the previous figures for 
a variety of values of d, the results of which are reported 
in Fig. [5 In particular, the long-time behavior of reff(t) 
after a quench at r = rc (such as the one displayed for 
d = 3.4 in the upper left panel of Fig. [5) can be fitted with 
the expected algebraic law in order to extract the 
value of a as a function of d. The resulting numerical es¬ 
timates are indicated by the dots in the upper left panel 




















of the figure, where they are compared with the theo¬ 
retical prediction in Eq. (411 (dashed line). While the 


agreement between the latter and the numerical data is 
very good for d < 3.6, slight deviations appear upon ap¬ 
proaching the upper critical dimensionality dc = 4 of the 
model, due to the expected corrections to scaling which 
are known to become increasingly relevant as d —>■ dc |63| . 

Analogously, by fitting the time dependence of the crit¬ 
ical iGxik = 0,t,t) with the algebraic law ^ G, one 
can estimate the numerical value of 7 as a function of 
d. These estimates, reported in the upper right panel of 
Fig. [^(symbols), can then be compared with the analyt¬ 
ical prediction 7 = 1 -f 2a = d/2 for d < 4 and 7 = 2 
for d > 4 (dashed line) which follows from Eqs. (311, 
(411, and ([4^. Also in this case the agreement between 


the data and the analytical prediction is very good apart 
from a region around dc, where corrections to scaling 
make the extraction of the exponent from the data more 
difficult. While both of the previous evidences in favor 
of the theoretical predictions of Sec. |III| are based on 
the scaling properties of Gk and of refi(t), an indepen¬ 
dent and more stringent test is provided by the analysis 
of the response function Gn, whose analytical form in 


Eqs. (271 and (291 does not involve unknown parame¬ 


ters (such as those which fix, instead, the amplitude of 
Gxik = in Eq. ([3l|). In particular, another esti¬ 

mate of a can be obtained by fitting \GR{k = 0,t,t')/t\ 
with C<i>Q,(t'/t) (see Eq. ( [4^ ) as predicted by Eq. (291. 
The resulting values of a are reported in the lower left 
panel of Fig. together with the analytical prediction 
(dashed line) of Eqs. ([dlj) and (42). Alternatively, one 


can estimate the value of the proportionality constant C 
by fitting \Gii{k = 0 ,t,t')/t| with G^ait'/t) where a is 
now fixed to the theoretically expected value reported in 
Eqs. ( |4l| ) and ( [4^ . The resulting numerical estimates 
are indicated by the symbols in the lower right panel of 
Fig. I together with the analytical prediction (dashed 
line) C = l/(2a) which follows from Eq. (29) and from 
the theoretical values of a of Eqs. (41) and (|42|. Also 


for the lower left panel, the agreement with theoretical 
predictions is good, except for values close to dc, while it 
is remarkably good for the lower right panel. 

Further below we argue that the exponent 7 = 1 -|- 2a 
which describes the algebraic behavior of (fc, t, t) both 
at short and long times (see Eqs. (31) and (32)) is the 
same as the one introduced in Ref. m in order to char¬ 
acterize the small-momentum behavior ~ of pkit) 
up to a cut-off k* ^ . The quantity Pk{t) corresponds 

to the average number of excitations with momentum k 
of the pre-quench Hamiltonian which are produced af¬ 
ter a so-called double quench, i.e., when the parameters 
of the post-quench Hamiltonian are restored suddenly to 
their initial values after a time t has elapsed from the 
first quench occurring at t = 0. The values of 7 which 
were numerically determined in Ref. |44| for a quench to 
the critical point in d = 3 and 4, i.e., 7 = 3/2 and 2 are 
in perfect agreement with our numerical estimates and 
analytical predictions reported in the upper right panel 


of Fig .0 


B. Quench below the critical point 

The Keldysh (Fig. upper panels) and the retarded 
(Fig[^ lower panels) Green’s functions Gk and Gr, re¬ 
spectively, for r > Tc are characterized by an oscillatory 
behavior, which denotes the presence of a finite length 
scale ^ in the model, set by the asymptotic value r* of 
reff(t). On the contrary, for r < Tc, the effective param¬ 
eter reff(t) turns out to decay to zero as ~ acot~^ (see 
the upper left panel of Fig. [^, i.e., with the same power 
law as at criticality r = Vc' correspondingly, Gr^k ex¬ 
hibit algebraic behaviors, which however differ from the 
critical ones. In fact, it is rather related to the phe¬ 
nomenon of coarsening which we discuss further below 
in Sec. |IV C[ As we discussed in Sec. |HI[ the value of 
Oco — as well as of a for r = Tc — can in principle be 
determined in such a way as to satisfy the self-consistent 
equations (36) and (39), which indeed admit two differ¬ 
ent solutions, reported in Eqs. (41) and ( [4^ . While the 
former correctly describes the observed behavior at crit¬ 
icality (see the evidence presented in Fig. [^, it is quite 
natural to expect the latter to describe the other possible 
scaling behavior, i.e., the one associated with coarsening. 
Numerical evidence of this fact is presented in Fig. In 
particular, the upper left panel shows the value of Cco 
(symbols) for various values of d, as inferred by fitting 
the corresponding numerical data of reg{t) for r < Vc 
with The dashed line corresponds to the the¬ 


oretical prediction reported in Eq. (45). Although the 
numerical data reported in Fig. prefer to a quench with 
r = — 3 < Tc, we have verified that these numerical es¬ 
timates are not affected by the choice of r < re- Note 
that while a at r = Tc as a function of the dimensionality 
d shows a marked change in behavior upon crossing the 
upper critical dimensionality dc (see the upper left panel 
of Fig. [^, being zero above it, this is not the case for 
Oco- Analogous consideration holds for the other quan¬ 
tities discussed further below, when compared with the 
corresponding ones at criticality. Heuristically this might 
be expected based on the fact that — as in the case of 
classical systems |65l IHH] — coarsening for r < rc is gen¬ 
erally driven by a different mechanism compared to the 
one controlling the behavior at r = rc, which is related to 
critical fluctuations and which is therefore affected upon 
crossing dc (see Sec. |IV C ). As a peculiar feature of Oco, 
we note that it vanishes for d = 3. 

As argued above and demonstrated by the curves in 
the upper left panel of Fig. Gxik = 0,t,t) grows al¬ 
gebraically both at r = Tc and for r < rc, in the latter 
case as ^ H™. The upper right panel of Fig. shows the 
estimates of jeo obtained by fitting the numerical data 
for G/f(fc = as a function of the dimensionality 

d. The dashed line corresponds to the theoretical pre¬ 
diction 7(.o = 1-1- 2aco with aco given by Eq. ( 45) . As 
it was done in Fig. for r = rc, the lower panels of 
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Fig. 0 consider G[i{k,t,t') at A: = 0. By assuming that 
the scaling behavior in Eq. (291 carries over to r < Tc, 
\Gii{k,t,t')/t\ is fitted by Cco^a^^{t'/t) (with given 
in Eq. ©) in order to extract Oco (lower left panel) or 
to estimate Geo once Oco has been fixed to its theoretical 
value in Eq. ( [4^ . In both panels the corresponding the¬ 
oretical predictions are reported as dashed lines and, as 
in the case of the upper panels, the agreement with the 
numerical data is excellent, with a hint of slight devia¬ 
tions upon approaching the lower critical dimensionality 
d = 2 of this model. 

As we mentioned at the end of Sec. |IV Aj the exponent 
7 co = 1 + 2aco discussed here in connection to the scaling 
of GK{k,t,t) (see Eqs. (311 and (32)) is the same as the 
exponent 7 introduced in Ref. |44| in order to characterize 
the scaling behavior of Pk{t)- The values of 7 which was 
numerically determined in Ref. |44| for a quench below 
the critical point in d = 3 and 4, i.e., 7 = 2 and 3 are 
in perfect agreement with our numerical estimates and 
analytical predictions reported in the upper right panel 
of Fig. 1^ 


C. Coarsening 

The numerical data presented in Sec. |IVB] clearly show 
that the non-equilibrium dynamics of the system after a 
quench to r < Tc features an emerging scaling behav¬ 
ior which we partly rationalized in Sec. |III| and which 
is characterized by scaling exponents depending on the 
spatial dimensionality d. As anticipated, these scaling 
forms are expected to be related to the coarsening dy¬ 
namics |421 H5] , analogously to what happens in classi¬ 
cal systems after a quench below the critical tempera¬ 
ture I15HS7] (see Appendix [A| for the discussion of a spe¬ 
cific example). In fact, when a classical system prepared 
in a disordered state is quenched below the critical tem¬ 
perature, the global symmetry cannot be dynamically 
broken and, consequently, the order parameter remains 
zero in average. Nevertheless, symmetry is broken lo¬ 
cally by the creation of domains within which the order 
parameter ip takes the value characterizing one of the 
possible different and competing phases. The average 
linear extension L{t) of the ordered domains increases 
with time t, until a specific domain possibly prevails over 
the others, establishing the equilibrium state. However, 
because of such competition. Lit) grows algebraically as 
L(t) oc where Zc > 0 is an exponent depending on 

the universal properties of the model, and equilibrium 
is reached only in an infinite time. Consequently, this 
lack of an intrinsic length scale in the system affects the 
equal-time two-point correlation function 

C{r,t) = {if{-x + r,t)ip{yi,t)), (49) 

and its spatial Fourier transform G{k,t), which, accord¬ 
ing to the scaling hypothesis [SS], are expected to display 
the scaling forms 

C{r,t) = f{r/L{t)) and G{k,t) = [L{t)fJ{kL{t)), (50) 


where d is the spatial dimensionality and f{x) the Fourier 
transform of f{x). 

The scaling forms for a quench below Tc highlighted 
in Sec. IVB do not satisfy the scaling hypothesis ( [50| , 
as it was noticed for d = 3 in Refs. |42l 1^ . In fact, 
the Keldysh Green’s function at equal times Gfcik,t,t) 
— which corresponds to the correlation function G{k, t) 
mentioned above — can be written as a scaling form by 
using Eqs. (34) and (|4^, which reads: 


GK{k,t,t) = [L{t)r^^gd{kL{t)), 


(51) 


where Gdix) is the scaling function, L{t) cx t (i.e., the 
coarsening exponent Zc takes the value Zc = 1) and jeo = 
d — 1. As this 7co differs from d, Eq. violates the 
scaling form (50) in all spatial dimensions. 


V. CONCLUSIONS AND PERSPECTIVES 

In this work we provided a complete characterization of 
the dynamical scaling which emerges after a deep quench 
of an isolated quantum vector model with 0{N) symme¬ 
try at or below the point of its dynamical phase tran¬ 
sition. The lack of intrinsic time and length scales is 
responsible for the occurrence of aging phenomena simi¬ 
lar to the ones observed in non-equilibrium classical sys¬ 
tems giino] or, more recently, in isolated m quantum 
many body systems. While previous investigations of 
this phenomenon were based on a perturbative, dimen¬ 
sional expansion around the upper critical dimensionality 
dc = 4 of the model @7], here we carry out our analy¬ 
sis within the exactly solvable (non-perturbative) limit 
IV —>■ 00, which allows us to obtain exact results for scal¬ 
ing exponents and scaling functions of the relevant dy¬ 
namical correlations, depending on the dimensionality d 
of the model. We find that the value of the pre-quench 
spatial correlation length (assumed to be small) controls 
the microscopic time Ia after which the aging behavior 
emerges: in addition, it acts as an effective temperature 
for the dynamics after the quench, which, inter alia, de¬ 
termines a shift of the upper critical dimensionality dc of 
the model, as it occurs in equilibrium quantum systems 
at finite temperature |681169| . Moreover, we provide evi¬ 
dence of the emergence of a dynamic scaling behavior for 
quenches below Tc, associated with coarsening, which we 
characterized numerically and analytically by studying 
the dependence of the relevant exponents on the spatial 
dimensionality d of the system. 

The exactly solvable model considered here provides 
a prototypical example of a dynamical phase transition 
(DPT) and of the associated aging occurring in a non- 
thermal stationary state. This state is expected to be¬ 
come unstable in systems with finite N, when the non- 
integrable terms of the Hamiltonian become relevant, 
causing thermalization. Nonetheless, this DPT might 
be still realized in the prethermal stage of the relax¬ 
ation of actual quantum systems evolving in isolation 
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from the surrounding environment [HllIMni- The lat¬ 
ter are nowadays rather easily realized in trapped ul¬ 
tracold atoms, the behavior of which can be analyzed 
with remarkable spatial and temporal resolution [Mia. 
In general, physical systems which can be described by 
some effective Hamiltonian with 0{N) symmetry include 
experimental realizations with ultra-cold atoms of the 
Bose-Hubbard model O [TJ IE] (corresponding to = 2) 
and one-dimensional tunnel-coupled condensates nana 
(N =1). In passing, we mention that an alternative 
and promising experimental realization of these models 
currently under investigation m is based on fluids of 
light propagating in non-linear optical media, which is 
expected to be ready for testing in the near future. 

At least in principle, the universal dynamic scaling be¬ 
havior emerging after a sudden quench which is high¬ 
lighted in the present work can be experimentally studied 
by determining directly the two-time linear response and 
correlation functions of the system. Alternatively, one 
can exploit the statistics of excitations produced after a 
(double) quench, as proposed in Ref. [H]. In fact, the n- 
th cumulant Cn (t) of the corresponding distribution was 
shown to grow as a function of the time t elapsed since 
the quench, with a behavior which may saturate, grow 
logarithmically or algebraically, depending on n, d and 
on whether the quench occurs above, at, or below criti¬ 
cality. In particular, it was shown |33] that the increase 
in time of Cn{t) is proportional to the integral over k of 
the n-th power of the quantity pk (t) related to the num¬ 
ber of excitations, which we briefly discussed in the last 
paragraph of Sec. |IV Aj In turn, at long times, the leading 
growth of pk{t) is the same as the one of |/k(t)P (see the 
definition of pk{t) in Ref. (B]), i.e., of iG (see 
Eq. (131 here). As a result, a simple comparison with 
Eqs. (311 and (l32|) yields 


C7„(t)(x y’d^fc[iG;f(fc,t,t)]” (52) 

for a quench to the critical point; for a quench below it, 
instead, one finds the same expression with a replaced 
by aco, i.e.. 


Cn(t) cx yd‘^/c[*Gif(fcG,i)]” 


(53) 


where the values of the exponents a and aco are given in 
Eqs. (41), (42), and (45). As a result, a measure of the 


statistics of the number of excitations for a quench would 
provide direct information on the aging and coarsening 
properties of the system. 

The quantum aging and coarsening discussed in this 
work enrich the list of mechanisms underlying the 
scale-invariant non-thermal fixed points (NTFP) |71H73| . 
which have been so far interpreted in terms of quan¬ 
tum turbulence IBHZS] and dynamics of topological de¬ 
fects diiis]. The extent to which these mechanisms are 
interconnected and combined in the dynamics of physical 
systems represents an intriguing yet challenging question 
for future investigations. 
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Appendix A: Comparison with aging after a classical 
quench 

In this Appendix, we briefly review the phenomena of 
aging and coarsening in classical systems evolving in con¬ 
tact with a thermal bath after quenching its temperature 
either at or below a critical point of the system (48]. In 
both cases, the lack of intrinsic time- and length-scales 
causes the emergence of algebraic behaviors in the tempo¬ 
ral dependence of, e.g., two-time correlation and response 
functions which can be characterized in terms of scaling 
exponents and scaling functions with a certain degree 
of universality. The stochastic dynamics of these classi¬ 
cal statistical systems can be simply described by effec¬ 
tive models m, which take the form of Langevin equa¬ 
tions for the (coarse-grained) relevant degrees of freedom 
of the system. For example, a A-component real held 
if = (ifi,..., pn) obeying purely dissipative dynamics 
evolves according to the so-called model A 


dt<Pai.^,t) = -D 


5'H[p\ 

Spa{^,t) 


+ Ca(x,t), (Al) 


where D is a diffusion coefficient, H a 0(A)-symmetric 
effective Hamiltonian in d spatial dimensions 


U[ip]= / d‘^x 


{Vpf 


.V’ 


4!A 




(A2) 


and Ca is a zero-mean Gaussian white noise describing the 
thermal fluctuations of the reservoir, with correlations 

(Ca(x,<)C6(x',t')) =2Z?T(5,b(5(‘')(x-x')5(t-t')- (A3) 


The evolution prescribed by Eq. ( |A1[ ) is such that the dis¬ 
tribution of the fluctuating field p at long times relaxes to 
the equilibrium distribution Peqlv’] oc indepen¬ 

dently of the initial condition. In this stationary state, 
and upon varying r, the system undergoes a second-order 
phase transition at the critical point r = Vc oi H. Gener- 
ically, r in the effective Hamiltonian (A2) is actually a 
function of the temperature T, with r — oc T — T^,, 
where is the critical temperature of the classical sys¬ 
tem. 

A classical quench protocol |4^ [50] consists in prepar¬ 
ing the system in, e.g., a disordered state at high tem¬ 
perature (i.e., with vanishing correlation length) at t = 0 
and in letting it evolve with fixed r = Tc or r < rc for 
t > 0. As a result, in both cases, the system relaxes 
to the equilibrium distribution with an algebraic behav¬ 
ior characterized by universal exponents which control. 
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e.g., the scaling of the Fourier transform in space of the 
two-time and two-point correlation and response func¬ 
tion C{k,t,t') and R{k,t,t'), respectively. In particular, 
for a quench to the critical point, C{k, t, t') and R{k, t, t') 
can be calculated exactly in the limit N ^ oo |491 [HO] 
and they read (by rescaling time one can set D = 1) for 
2 < d < 4: 

/ . \ (4-d)/4 

R{k,t,t')=e{t-t')^^^ (A4) 

and 

C{k, t, t') = ^F{k^t, (A5) 

with F[x^y) a scaling function defined as F{x^y) = 

Assuming for 

simplicity t > t', the asymptotic expression of C{k,t,f) 
at short times t, t' ^ k~^ is 

4 f t\ 

, (A6) 

while, at long times t, t' ^ k~‘^, it reads 


correspond, respectively, to the classical R and C — turn 
out to agree with Eqs. (A4|, (A6l, and (A7). 

The effect of having an isolated instead of an open 
quantum system is not only revealed by the different 
value of the dynamical exponent z (1 and 2, respec¬ 
tively) but also by the scaling form of the correspond¬ 
ing correlation and response functions at short times, 
where the actual value of z does not appear explicitly. 
In fact, from Eqs. p4| , ra, and p0| one finds that 
iGK{k = ~ which can be compared with 

C{k = 0,t,t') in Eq. ( |A6| ). While the dependence of the 
earliest time t' is characterized by the same exponent d/4 
this does not apply to the dependence on t and in fact 
the overall scaling form is significantly different, as ex¬ 
pected on the basis of the different scaling dimensions 
of the relevant fields IIZIIIS]. Comparing, instead, the 
corresponding short-time response functions in Eqs. (291 
(see also Eq. ( [4l] )) and (A4|, one finds that the depen¬ 
dence on the time f is characterized by opposite powers, 
whereas the ones on t are seemingly unrelated. 

The classical model can be exactly solved in the limit 
N ^ oo also if the quench occurs from the disordered 
state to below the critical point r < [50]; for d > 2, 

the corresponding response function reads: 




(4-d)/4 




(A7) 


The effective classical dynamics prescribed by Eq. (Al) 


is actually relevant also for the quantum system with 
Hamiltonian 0 investigated in the previous sections, 
when the field cf) is linearly coupled to a bath of harmonic 
oscillators with ohmic spectral density and at equilib¬ 
rium with temperature T. This open quantum system 
is then characterized by an equilibrium critical point at 
r = rl'^{T) at which the microscopically large-distance, 
long-time properties are effectively described by the criti¬ 
cal classical model described above |58|. In fact, a critical 
quench in this open quantum system was recently stud¬ 
ied in Ref. starting from a disordered pre-quench 

state and the scaling functions of Gr and Gk — which 


t \ 


d/4 




R{k,t,t') = e{t-t') j 

while the correlation function is (with t > t') 

d/4 


(AS) 


C{k,tX) = 


(/) 




(A9) 


with Mgq being the value that the order parameter would 
have in equilibrium in the system described by the post- 
quench classical Hamiltonian. Accordingly, this system 
exhibits coarsening in the strict sense |65|, as Eq. ( |A9| ) 
satisfies the dynamical scaling in the form indicated by 
Eq. (501 with the proper value Zc = 2 of the coarsening 


exponent. 
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FIG. 1. (Color online). Effective parameter red as a fnnction of time t. Upper left panel: res for a quench above (uppermost 
curve), at (lowermost curve), and below (intermediate curve) the critical point rc with d = 3.4, = 10^, u ~ 80, and a 

Gaussian cut-off function h{k/A) = exp[—fc^/(2A^)] with A = tv/2. The black dashed lines indicate an algebraic decay ~ 
Upper central panel: r^fi at r = Tc for sharp (uppermost curve) and smooth (lowermost curves) cut-off functions h{k/A) 
indicated in the legend (see also the main text) with d = 3, SIq = 10^! w ~ 47, and A = 7r/2. Upper right panel: reS at r = rc 
for various values of the sharp cut-off A,increasing from top to bottom, with d = 3, Hq = 10^, and u ~ 4.74, compared with 
the expected algebraic decay ~ (dashed line). Lower left panel: refr at r = rc for various values of the initial parameter floi 
decreasing from top to bottom, with d = 3, u ~ 4.74, and a Gaussian cut-off function h with A = 7r/2. Lower central panel: 
reff at r = rc for several values of d > dc = 4, w = 1, flo = 10^ and a Gaussian cut-off function h with A = 7r/2; the dashed 
lines superimposed on the curves are proportional to (uppermost curve), (intermediate curve) and ® (lowermost 

curve). Lower right panel: reff at r = rc for several values of flo,decreasing from top to bottom, with fixed d = 4.5, u = 1 and 
a Gaussian cut-off function h with A = 7r/2. 
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FIG. 2. (Color online). Keldysh and retarded Green’s functions Gk and Gr, respectively. Upper left panel: GK{k,t,t) at 
k = 0 and equal times as a function of time t, for r below (uppermost curve), at (intermediate curve), and above (lowermost 
curve) the critical value Tc. The dashed lines superimposed to the curves for r = and r < are proportional to and 
, respectively. Upper right panel: GK{k,t,t) at fixed equal times t = 10® as a function of momentum k after a quench 
with r below (yellow, uppermost line), at (red, intermediate line) and above (blue, lowermost line) the critical value Vc- The 
dashed lines superimposed on the curves for r = Tc and r < Vc are proportional to kT^^'^ and fc“®, respectively. Lower left 
panel, main plot: retarded Green’s function \GR{k,t,t')/t\ for k = 0 and t = 2 x 10® as a function of t'/t for quenches of r 
at (lower curve) and below (upper curve) the critical value r^. The black dashed lines superimposed on the curves for r = rc 
and r < Tc correspond to the theoretical predictions 2$i/4(t'/t) (see Eq. (47l) and /t), respectively. Inset: retarded 

Green’s function GR{k, t, t') for k — 0 and t — 2x 10® as a function of t' for quenches of r above the critical value r^- The black 
dashed line corresponds to the prediction in Eq. (|30| with k —>■ y/r* (see the main text) and r* ~ 0.34. Lower right panel: 


retarded Green’s function GR{k,t,t') as a function of k{t — t') for t = 2 x 10®, t' = 1.9 x 10® and r below (main plot, yellow, 
superimposed on the red one), at (main plot, red, superimposed on the yellow one) and above (inset, blue, solid curve) the 
critical value r^. In the main plot, the two curves are indistinguishable from the superimposed dashed line, which corresponds 
to the theoretical prediction in Eq. (30l. In the inset, the curve is indistinguishable from the superimposed dashed line which 
corresponds to Eqs. (301 and (481 with r* ~ 0.34. In all the panels of this figure, the spatial dimensionality is fixed to d = 3, 
= 10®, u ~ 47, and a Gaussian cut-off is used with A = 7r/2. 
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FIG. 3. (Color online). Nnmerical values (symbols) of the exponents 7 , and a and prefactors a and C as functions of the spatial 
dimensionality d for a quench at the critical point. Upper left panel: coefficient a of the effective parameter reff(t) computed 
numerically compared with the theoretical value a = d{l — d/‘i)l4: (dashed line). Upper right panel: exponent 7 obtained by 
fitting iGKik = 0,t, t) with an algebraic law oc G■, the dashed line indicates the analytical value 7 = d/2 for d < 4 and 7 = 2 
for d > 4. Lower left panel: exponent a obtained by fitting |Gij(fc = 0,t = 2 x lQ^,t')/t\ with /t), see Eq. (47l; the 

dashed line represents the theoretical value a = (d — 2)/4 for d < 4 and a = 1/2 for d > 4. Lower right panel: prefactor C of 
GR{k — as a function of the dimension d, obtained by fitting \Gi{{k = 0,t = 2 x 10^,t')/t\ with G^a{t'/t), with a given 

by the theoretical values rep orted above; the dashed line represents the theoretical prediction for G = 2/(d — 2) for d < 4 and 
G = 1 for d > 4 (see Eq. (291). The numerical data presented in this figure have been obtained with a Gaussian cut-off function 
with A = 7r/2 and Hq = lO*, while for each point a different (inconsequential) value of u has been used. The statistical error 
bars on the numerical points are smaller than the symbol size. 
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FIG. 4. (Color online). Numerical values (symbols) of exponents 7co, and a^o and prefactors a^o and Geo as functions 
of dimension d for a quench below the critical point. Upper left panel: numerical estimates for the coefficient a^o of the 
algebraic decay of the effective parameter res{t) — as a function of dimension d, compared with the theoretical value 

dco = (3 — d)(d — l)/4 (dashed line) . Upper right panel: numerical estimates of the exponent 7co obtained by fitting 
iGif(fe = 0,t,t) with the dashed line represents the analytical prediction 7co = d — 1. Lower left panel: exponent Oco 
obtained by fitting \Gii{k = 0,t = 2 x 10^,t')/t\ as Geo “haco (^V^) (with $a(a;) given in Eq. @); the dashed line represents the 
theoretical value aco = {d— 2)/2. Lower right panel: prefactor Geo of GR{k,t,t') as a function of the dimension d, obtained by 
fitting |Gfl(fc = 0,f = 2 X lQ^,t')/t\ as Gco^a^o{t'/t), with a^o fixed to its corresponding theoretical value (see Eq. (45l); the 
dashed line represents the theoretical value Geo = l/(d— 2). The numerical data presented in this figure refer to a system with 
flo = 10"^! ^ = ~3, and a Gaussian cut-off function with A = 7r/2; for each point a different (inconsequential) value of u has 
been used. The statistical error bars on the numerical points are smaller than the symbol size. 







